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SUMMARY 

A new computer code has been developed to analyze chemically reactive 
flow and spray combustion processes occurring inside a strati f ro y 
eng'nr Mathematical and numerical details of the new code ■« « " 
described by the present authors. This paper presents the results of limited, 
initial computational trials as a first step in a long-term assessment/ 
validation process. The engine configuration studied was chosen to a PP^^ [ I iate 
existing rotary engine flow-visualization and hot-firing test rigs. yp 
results include (1) pressure and temperature histories, (2) torque genera 
bv the nonuniform pressure distribution within the chamber, (3) energy please 
rLpf and 1 41 various flow-related phenomena. These are discussed and com- 

pared' with ither predictions reported in the literature^ ee J h f 0 fi Corporal inf 
for improvement in the spray/combust ion models and the need for incorporat g 

an appropriate turbulence model are also discussed. 


INTRODUCTION 

The rotary combustion engine (RCE) would be a des ! rabl ®. p "^P‘ a "‘ d f ^ dus _ 
light aircraft, drones, auxiliary and ground power units, and marine and Indus 
trial aDDlications if its efficiency could be made closer to that of diesel 
engines. P It has inherent advantages over reciprocating engines such as ^ ghe 
aif no. capacity, higher power-to-.eight ratio, ln«-n aa reaandles 
vibration. An initial attempt to introduce a gasoline-fueled rotary engine 
into the general aviation market in the mid 1970 s was unsuccessful , however, 
because of poor fuel economy, uncertain aval labi lity of avgas and * r<jh 

weight advantage over contemporary reciprocating engines. ^ . <• 

sponsored by industry, NASA, and the Navy has led to the development of the 
straUf ied-charge roUry engine (SCRE) concept. This work has a >- a dy demon- 
strated a substantial improvement in specific P° wa r , mul t ifue capability, ana 
a reduction of cruise brake specific fuel consumption (BSFC) from over 

0.5 lb/bhp-hr to a value approaching 0 •f 1 P ! b/ $ h RT hr ' • A Methods to 

improvement was the use of Computational Fluid Dynamics (CFD ) methods 
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analyze the airflow, spray, and combustion events nf thp rpf • 

«uise C BSFC d from'lh" 6 " 1 Sp °f 0r ? d by N4SA is “i""*) at further reducing^ 
cruise B5FC, from the current value to 0.35 or less bv thp pnH nf iqqo u u 

tion snrp eCte H im ^ vement " iIJ enabled by further CFD-drlJen fj,el ini ec 

and related changes. Thus, the ava lab ? t^f lcc“a?e feHab e 
opera i ingCFD simulations of the RCE internal processes is c ^v’a fev e l 
in the future development of this promising eng?™ ?oncep? Y 

models a Uefs° d ?‘and " ank 5! e " 8 l" e * ere ba ^ d on thermodynamic 

SSS'XI"' teL,*!* 'a?!" 

2% “!ytX C ?S? e ?„!e , ,he r °’? ry en8ine at“ee?i J3 &.T 

U } M le? e K,V ? COde deV? * 0Ped ai Ss'M^:' N aU™aI F Labo 8 ?a,;? y a 

»S£C- 

Trpf ft i nro a l ime w ^ en en gme speed becomes small. Shih et al 
Engine n thrabsir^f irS R t T dimen Ti° nal com P ut ations of a motored Wankei 
Bea. .arming type of al.ernating-diren^ri^i^if^J^^L^r'cLpu 

(ref 9) Linea?'1tabUitv X ln2| ed - ,0 h ,hre n dime C sions »" Steinthorsson et al. 

tionallv stah?p ?n !wn H- ty y! V S has shown th »t the ADI method is uncondi- 
y stable in two dimensions but is uncondi tional lv nnstahip in 

elce«!?? S ' 41 : h ou g h artificial dissipation has s«e s?ab??,z?i e?f' c [ a „ 

Recently L?°et ty an d reduce accuracy and convergence. 

schemes ioMth^ iith ^ 1 " ,od,f : ed ,h ?' r I-EWIS-3D code based on up»ind 

together *ith the incorporation of a k-e turbulence model. 

present autho^trpffi L h C ° de were recentl V described by the 
Hrnfe^^? ^ ^ aS *°^ 8aUra ^ aP ^ eaaa ^ a ^ Vaa * eaa °*® n_ ^ na ^®tms*of ^ressur^and 

h..t?nn atU ?a hl fi! 0ri k S V l0rqUe 8 enerated due to the nonuniform pressure distri 
bution •! thin the chamber, energy release rates, and various no«???Uted nhT 

in the iite?ati?e USS The'nh'“”? ariSOn ? Hh °! her P redicti o"a reported to da?e 

improvement in the spray/cimbus??o? S mode?rand , ?o t a?se?? q ?h? y „?e r d n ?o? f^corno 
rating an appropriate turbulence model. d f ° r lncor P°- 


PHYSICAL DESCRIPTION 


A schematic of the Wankei engine is shown in figure 1 The Wankei eneinp 
18 C ° mpOSed ° f 3 trochoid housing (ref. 1 ) with perigeral inlSke ln6 ^exhau 8 " 
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ports, fuel injector and spark igniter, a three-flank rotor, and a crankshaft. 
The contour of the inner surface of the outer casing of the Wankel engine is a 
two-lobe peritrochoid. The contour of a rotor revolving along an outer housing 
is a peritrochoid inner envelope. The rotor surface is modified to include 
rotor pockets (refs. 2 to 4). The presence of the rotor pockets not only 
alters the expansion and compression ratio of the engine, but also plays an 
important role in modifying its internal flowfield, mixing, and combustion 
characteristics . 


The present rotor configuration is adopted from Steinthorsson et al. 
fref. 9). In its other dimensions, the engine has a generating radius, R, ol 
0.1064 m; eccentricity, E, of 0.01542 m; clearance, C, of 0.004 m; chamber 
width, W, of 0.07 ra; and port width, W p , of 0.05 m. These dimensions give rise 
to maximum and minimum volumes of 750 and 115 cm 3 , respectively, thereby yield- 
ing a compression ratio of about 6.5. 


The rotor turns eccentrically at one third of the crankshaft speed. The 
three combustion chambers of the Wankel engine are the three regions enclosed 
between the three rotor faces and the peritrochoid housing, two side housings, 
two side seals, and lead and lag apex seals. In the present calculations, only 
one of the three combustion chambers is considered, since leakage through the 
seals is assumed to be negligible. As the rotor revolves around the crank- 
shaft, each of the combustion chambers is continually deformed. This produces 
the necessary compression and expansion of the fluid for the required engine 
performance during each cyclic operation. 


It should be noted that each combustion chamber has two minimum and maxi- 
mum volume positions. Hence, the top-dead-center (TDC) and bottom-dead-center 
(BDC) positions are defined differently than would be the case in a reciprocal, 
ing engine, and it will be important to keep the following conventions in mind. 
In this study, a crank angle of zero radians corresponds with the minimum vol- 
ume that occurs near the ports. This position is referred to as the BDC. e 
other position, corresponding with the minimum volume near the tuel injec or 
location, is referred to as the TDC. Thus, combustion occurs near TDC as in a 
reciprocating engine, but the BDC position has a different meaning. 


The exhaust port opens at a crank angle of 5.96 rad before BDC and closes 
at 1.07 rad after BDC, and the exhaust back-pressure is maintained at 0.85 atm. 
The intake port opens at 1.26 rad before BDC and closes at 5.96 rad after BDC. 
Fresh air enters through the intake while the intake stagnation temperature 
Tc int and static pressure Pint are maintained at 400 K and 1.25 atm, 
respectively. The fuel injector has eight holes and is located near the middle 
of the trochoid housing as shown in figure 1. The fuel emerges in a fan-shaped 
pattern consisting of eight sprays. The initial spray distribution is given 
bv a Rosen-Ramler distribution with the initial droplet sizes ranging trom 
10 um < rk i < 30 pm. The initial droplet velocity and temperature are taken 
to be 100 m/s and 300 K, respectively. The fuel injection begins at 8.6 rad 
after BDC and ends at 9.25 rad after BDC. The liquid fuel injection is simu- 
lated by injecting a specified number of discretized parcels of fixed mass 
associated with a given droplet-size group at uniform intervals corresponding 
with the fuel-injection time step. The ignition process begins at 8.65 rad 
after BDC and ends at 8.75 rad after BDC. The ignition spark is simulated by 
raising the temperature of a few nodes of the computational domain on the wail 
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at the spark ignitor location from 400 to 1500 K. A complete physical descrip 
tion of the engine is provided in our earlier paper (ref. 11). 


RESULTS AND DISCUSSION 

Numerical solutions were obtained for the engine flowfield under motoring 
conditions for a single case with all of the ports closed; they were also 
obtained under firing conditions for two cases, with equivalence ratios of 
0.45 and 0.75. For all cases, the computations were initiated before the open- 
ing of the exhaust port for only one of the three chambers (formed with the 
second rotor flank as shown in fig. 1). Here we present the results obtained 
for the three different cases for an engine speed of 4000 rpm. The tempera- 
tures of the rotor and housing surfaces were maintained at 300 K for the case 
considered under motoring conditions and at 400 K for the firing conditions. 

The initial conditions correspond to the conditions of quiescent air at a pres- 
sure, Pj n , of 1 atm and a temperature, Tj n , of 300 K for the motoring engine 
and at 1 atm and 400 K for the firing engine. 

All of the computations were performed with a uniform time-step size such 
that 15 000 time-steps span the 6t r radians of one entire rotor cyclic— rotation. 
The calculations were performed on a grid with a mesh size of i = 31, j = 16, 
and k = 20, where i, j, and k represent the coordinate surfaces in the 
direction extending from the trai ling-edge surface to the leading-edge surface 
of the combustor, from the rotor to the housing surface, and from the side wall 
to the symmetry plane of the domain between the end-to-end side walls, 
respectively. 


Motoring Results 

The usefulness of any numerical scheme lies in its ability to generate 
meaningful predictions for a wide variety of engine operating conditions. Both 
the numerics and the models used for turbulence, combustion, and sprays deter- 
mine the accuracy of the predictions. The selection of a finite-difference 
formulation is mainly dictated by a careful balance between the numerical sta- 
bility, efficiency, and accuracy considerations for the problem under consider- 
ation. Stability is often achieved at the expense of accuracy and efficiency. 
It is also instructive to note that a numerical scheme is often made stable 
either by explicitly adding some amount of artificial dissipation, or by the 
numerical dissipation that inherently arises from a given finite difference 
formulation. For the predictions to be meaningful, the numerical dissipation 
required for stability considerations should be small compared with the dissi- 
pation arising from the actual physical processes governing the flowfield. 

The best way to assess the accuracy of a numerical scheme is to compare 
the predictions with experimental data. In our previous paper (ref. 11) a good 
measure of qualitative agreement was noted between the predictions and an 
experimental flow-visualization pattern during the intake process. In the 
absence of applicable experimental data on the average pressure and temperature 
variation versus crank angle, a numerical experiment was performed for a motor- 
ing engine with all of the ports closed. This allowed comparisons between the 
predictions and the corresponding isentropic results to be made. 
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Figure 2 sho.sth. Mjyutjd --- r ra 8 ed e «e T pera , ure and volume^averaged 
resul t s^b t a i ned° from “Intropic cond iUons » -Ponan, of 

by the following reasons: 

this position approach 50 m/sec. 

(2) Viscous effects become increasingly significant near ™ s i nce„ 

surface area to volume ratio becomes ^ery arge resulting from the non- 

iStSTp^dii nb Uon Se &\ SlfffffS a unique feature of the 
rotary eSgfnls that has no counterpart in reciprocating engines. 

(3 ) The heat losses to the .alls also contribute to the observed differ- 
ences since the temperature of the nails is maintained at 300 K. 

m Because of the shape of the combustion chamber, all the gas within it 
may nit undergo the same degree of expansion or l““uk“i to CideJgo 

i"piiiion"ii l *hi 0 rotii^™eriroind°TDC This process enhances the convective 
heat transfer to the walls. 

,5, Numerical dissipation^ also ‘ t^con tribute^. the observed^ 
numer i ca^schemes 6 studied for a two-dimensional "“firing engine - ^ave not 

attempted to make any direct comparisons with their predctons be 

rsMis r fe.sj'.rr ^ 

about 3 percent. Even though our governing equat.ons^re^o condi _ 

strong-conservation law ( : . .ildj fference formulation weakly conserva- 

iive S a Si.i:eii”tii1bs^rod discrepancy in mass could be reduced further by 
increasing the number of grid points in the computational domain. 

Figure 3 shows the computed torque about the cen Jo a the J“'°' s l ”ner- 
K id 'insulting torque 

SZ -"«! C i-ii S‘iieds to b! performed -o-rcome^U torque ( f ^ 
iiS that h the°pressure Si^Tmoi "uniform near TDC and in fact . the 

'?T ^rq^edTii^ui -w» 

and -0.1259 Nm (seals) . 


5 



flnwfiplH °thf e n a bet er understanding of the observed nonuniformities in the 
it + ’ tbe P re f sur e and temperature contours together with the flow veloc- 

ity at the symmetric plane k = 19 and crank angle of 19.3 rad are shown in 

IhoJTas'do ?ed n^TH° Ur Pl ?.! S ’ - he isocontour lin es near the maximum are 
. otted lines and near the minimum are shown as solid lines. In the 

Hi cr* ty V h G I? r P °-’ only thre e different sizes of arrow symbols are used to 
* ab 1 t »* e variatl0n between the maximum and minimum values. Both the 

ajerseai an and e thp r fln re ? maximum vaIu f in the region near the trailing 

movement ^ ' d h f Vel ° Clty 15 essentially determined by the rotor 


Firing Results 


, * n 4 i i8 f eCtion we Present the results for two other cases that were con- 
ducted under firing conditions, for which the operating parameters were 
described earlier in the paper. The fuel/air equivalence ratio? 2, X taken to 
be 0.45 for one case and 0.75 for the other. 

Here we first present the results of the three-dimensional computations 
escribing the overall behavior of the combustor in terms of the volume- 

?a?iation^f S m?s? VGrSUS S°T buSt0r volume (P~ v diagrams). This includes the 

i tv lid v ?nnH?f Hn 8 S d te J" perature ’ the tor que due to pressure nonuni form- 
ity, and vaporization and combustion rates versus crank angle Then we D ro- 

romh..st eXamine tHe fl °Ti ie i d in d6taiI during the fuel vaporization and ? 
combustion processes. The details of the flowfield during the opening of the 

exhaust and/or intake are presented in detail in our earlier paper (ref. 11). 

ini .p? e computatio u s are initiated just before the exhaust port opens. The 

respei wi?? SUr Af?p? araper ? tur f of the § as ar e taken to be 1 atm and 400 K, 
respectively After the exhaust port opens the residual gas moves out of the 

combustion chamber, since the imposed pressure of 0.85 atm in the exhaust 

iwli nS Bith r thJ an the interior eng * ne Pressure during most of the compression 
[ Boh the en S ine Pressure and temperature remain fairly uniform during 
the exhaust process, as shown in figure 5. As the rotor moves farther the 

intalf 6 P ° rt opens anc * *-resh air moves into the combustion chamber. Most of the 
IV.I n P Z 6SS occurs during the expansion stroke of the engine. The slight 

pressure during the intake process merely reflects the intake § pres- 
sure, which is maintained at 1.25 atm. p es 

sion lltil rimhnc! ake u° rt closas \ the £ as undergoes nearly isentropic compres- 
sion until combustion begins. During the combustion process, both the pressure 

vah,es"near TDC "n* ? '° ab ° U ' 33 "i a " d 1450 K ’ »•**■« tW.’ZSS 

values near TDC. During the expansion stroke, the gas undergoes nearly isen- 

SS * he , exha “ S ! p ? rt ° pens a « ai " Both tS - InT 

redSceft^ „p f h 1 ^ * nterior Pressure near the exhaust port 

?hic pp a f below the imposed exhaust back-pressure of 0.85 atm. When 

h s condi ion is reached near a crank angle of 16 rad, the temperature remains 
p fa \l\l aaifo ™ at 650 K a brief period before the intake charge at 40 Q K 

thl nn the ch f lbar -. Th f temperature falls to 400 K during the remainder of 
the opening of the intake port. 

The torque variation versus crank angle in figure 6(a) shows that the 
pressure distribution on the rotor surface undergoes a kind of beating phenome- 
non associated with acoustic resonance during the intake and/or exhaust ? por t 
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opening period. It is also noteworthy that positive torque is generated during 
the overlapping period in which both the exhaust and intake ports are open 
simultaneously. Once again, the maximum negative torque is generated near the 
minimum volume as in the motoring case. The corresponding values ot the 
torque, averaged over one cycle, are found to be -0.817 Nm (rotor) and 
0.508 Nm (seals). Comparison with the previous case shows that the average 
rotor torque due to the pressure nonuniformity under both motoring and firing 
cases has a similar value. This indicates that the resulting pressure gradient 
on the rotor surface is a result of the fluid motion associated with the rotor 
movement. However, the average torque due to pressure acting on the seals tor 
a firing engine is found to be positive whereas it had a negative value tor the 

motoring case. 

Figure 6(b) shows the crank angle variation of the evaporated fuel and 
also the amount of reacted fuel. The total amount of liquid fuel injected is 
equal to 0.23 g, which corresponds to an equivalence ratio of about 0.45. lne 
slope of these curves indicates that most of the fuel, after it evaporates, 
reacts quickly with oxidizer to form products. This in turn implies that most 
of the fuel burns as if in a premixed-f lame environment. Although the validity 
of the results can be verified only after proper experimental data is made 
available, the extremely fast combustion rates observed in the present study 
are in contradiction to the numerical predictions presented by Abraham et at. 
(ref. 5), in which they reported lower rates for both fuel vaporization and 
combustion. Our present computations are based on the assumption of laminar 
fluid motion and chemical kinetics. In a turbulent reacting flow, ^e combus- 
tion rates are normally rate-controlled by the mixing rate rather than the 
Arrhenius reaction rate (ref. 12). Also one expects to predict lower tempera- 
tures for the gas near the walls, when the effect of turbulence is properly 
accounted for in the heat transfer calculations. This, in turn, might reduce 
both the vaporization and combustion rates. 

The effect of increasing the equivalence ratio from 0.45 to 0.75 is shown 
in figures 7 and 8. Although the qualitative behavior for both cases remains 
the same, the increase in heat input, as a result of increasing the fuel con- 
tent, is shown to result in an increase in the maximum temperature and pres- 
sure observed from 1450 to 1750 K and 33 to 44 atm, respectively. 


Flowfield Phenomena 

We now turn our attention to the details of the flowfield at the beginning 
of, during, and after the combustion and vaporization for the case with an 
equivalence ratio of 0.75. Figure 9 shows the isocontours of pressure and tem- 
perature, and the velocity vector plot at a crank angle of 8.6 rad. 1S I® 
before the fuel injection begins. Both the pressure and temperature distribu- 
tions remain fairly uniform in the region extending all the way from the trail 
ing apex seal to the end of the rotor pocket near the leading apex seal. The 
maximum values for the pressure and temperature are observed to occur within 
this region. Fairly large gradients in the distribution of both temperature 
and pressure are also observed in the clearance region near the leading apex 
seal. The fluid motion is shown to be strongly influenced by the rotor 
movement . 

Figures 10 to 12 show the droplet trajectories; the isocontour lines of 
fuel, oxidizer, carbon dioxide, and temperature; and the velocity vector plot 
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at a crank angle of 9.3. The figures showing the rotor pocket cross section 
are drawn in an elongated scale compared with the others in order to provide a 
better graphic illustration of the results. The polydisperse character of the 
spray is represented by different sized circles which are indicative of the 
size of the initial droplets. Although the particles retain the fan-shape con- 
figuration during this later part of the fuel vaporization process, the deflec- 
tion of the particles in the direction of the gaseous flow is evident from the 
path originally dictated by the particle initial conditions. 

A very small region near the fuel injector location is found to be fuel 
rich. Whereas stratified charge gives rise to a gaseous diffusion flame, a 
careful examination of figure 10 also reveals an absence of fuel concentration 
in the region near the rotor pocket surface where fuel concentration is 
expected otherwise from the presence of liquid fuel in that region. It is 
more likely that the combustion characteristics in that region are influenced 
by an isolated combusting droplet behavior. The oxygen concentration is found 
to be lower near the high temperature region than in the regions near the 
seals. An opposite trend is observed for the carbon dioxide distribution. 

The temperature in the flame region is around 3106 K and on the walls is 400 K. 
Again the fluid motion is found to be essentially dictated by the rotor motion. 

During the expansion stroke, after the combustion is completed, the dis- 
tributions of temperature and of oxidizer and carbon dioxide concentrations are 
found to become fairly uniform throughout the chamber as shown in figure 13. 


CONCLUDING REMARKS 

We have presented preliminary results for a Wankel engine under both 
motoring and firing conditions, obtained from the solution of unsteady, three- 
dimensional Navier-Stokes equations. These results are based on the assumption 
of laminar fluid motion and chemical kinetics, with the use of appropriate sub- 
models for combustion and sprays. 

The results show the average rotor torque generated by the nonuniform 
pressure distribution is nearly the same for the different cases considered 
under both motoring and firing conditions. This indicates that the pressure 
nonuniformity is determined mostly by the rotor-induced fluid motion. The 
pressure is higher near the trailing apex region than near the leading apex 
region . 

One indication of our study is that, under the present assumptions, vapor- 
ization appears to be more rate-controlling than mixing during the combustion 
process. This finding is contrary to the numerical predictions reported by 
Abraham et al. (ref. 5). While considering the fact that the present engine 
configuration is somewhat different from the one used in their studies, it 
should be noted again that the present computations are performed without 
attempting to describe the turbulent aspect of the fluid motion. Both the heat 
transfer, fluid mechanics, and combustion characteristics in a engine flowfield 
will be influenced strongly by the turbulent fluid motion. Our solution proce- 
dure is currently being modified to include an appropriate turbulence model. 

The code takes about 10 CPU hours on a CRAY Y-MP, when the calculations 
are performed on a 31 by 16 by 20 grid. This covers one entire cyclic period 
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of 6 it rads for a firing case. It appears that we are the first to have 
reported numerical predictions obtained from three-dimensional computations 
in terms of pressure, temperature, and torque histories for an entire cyclic 
period. The code appears to be potentially very efficient compared with the 
other schemes used in modeling Wankel engines. 
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(a) Pressure contours (P max = 0 77x10 6 , P min = 0 74x10 6 , and AP - 0.21x10 ). 



A Minimum 


(b) Temperature contours (T ma x = 479 4, T m j n - 300, and AT 5) 



(c) Velocity vector plot (V max - 43.8 m/s). 

Figure 4. — Gas pressure, temperature, and velocity at 0 = 19.3 rad and K = 19 with all ports closed 



35x10 s 



Volume, cm3 

(a) Pressure versus volume. 



(b) Temperature and volume versus crank angle. 

Figure 5.— Pressure, temperature, and volume histories forO = 0.45. 
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Volume, cm3 


Torque due to pressure on the rotor, Nm 



45x10 s 



Volume, cm3 

(a) Pressure versus volume. 



Figure 7. — Pressure, temperature, and volume histories for O = 0.75. 
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Volume, cm3 



Torque due to pressure on the rotor, Nm 



(a) Torque versus crank angle. 



Figure 8,— Torque and fueling histories with O = 0.75. 


Torque due to pressure on the seals, Nm 


Maximum 



(a) Pressure contours <P ma)( = 0.80x10 6 , P min = 0.76x10 6 , and AP = 0.26x10 4 ). 



(b) Temperature contours (T max = 606 1, T min = 400.0, and AT = 15). 



(c) Velocity vector plot (V max = 35.34 m/s). 

Figure 9 —Gas pressure, temperature, and velocity at 0 = 8.6 rad and K = 9 for the case with <r> = 0 75 
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C max = 0.0679, C min = 0, AC = 0.005 
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(a) Droplet trajectories. 




(b) Fuel mass fraction contours. 


Figure 10. — Droplet trajectories 


and fuel mass fraction contours at e = 9.3 rad for <I> - 0.75 










(a) Oxygen concentration contours (C max = 0.0923, C min = 0.0895, 
and AC = 0). 


(c) Temperature contours (T max = 2017, T min = 400, and AT = 116). 



(b)Carbon dioxide concentration (C max =0.1194, Cmin =0.1171, 
and AC = 0). 


(d) Velocity vector plot (V max = 21.85 m/s). 


Figure 13 — Oxygen, carbon dioxide, and temperature contours and gas velocity at 0= 11 rad and K = 19 for <t> = 0.75. 


22 



NASA 

National Aeronautics and 
Space Administration 

1. Report No. 

NASA TM- 103 104 


4. Title and Subtitle 


Report Documentation Page 


2. Government Accession No. 


Computational Experience With a Three-Dimensional Rotary 
Engine Combustion Model 


3. Recipient's Catalog No. 


5. Report Date 


6. Performing Organization Code 


7. Author(s) 

M.S. Raju and E.A. Willis 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


8. Performing Organization Report No. 

E-5249 

10. Work Unit No. 

505-62-11 


11. Contract or Grant No. 


13. Type of Report and Period Covered 
Technical Memorandum 


14. Sponsoring Agency Code 


15. Supplementary Notes 

Prepared for the Joint Symposium on General Aviation Systems cosponsored by the AIAA and FAA, Ocean City, 
New Jersey, April 11-12, 1990. M.S. Raju, Sverdrup Technology, Inc., Lewis Research Center Group, 

2001 Aerospace Parkway, Brook Park, Ohio 44142. (Work funded by NASA Contract NAS3-25266.) 

E.A. Willis, NASA Lewis Research Center. 

16. Abstract 

A new computer code has been developed to analyze the chemically reactive flow and spray combustion processes 
occurring inside a stratified-charge rotary engine. Mathematical and numerical details of the new code were 
recently described by the present authors. This paper presents the results of limited, initial computational trials as 
a first step in a long-term assessment/validation process. The engine configuration studied was chosen to approxi- 
mate existing rotary engine flow-visualization and hot-firing test rigs. Typical results include (1) pressure and 
temperature histories, (2) torque generated by the nonuniform pressure distribution within the chamber, (3) energy 
release rates, and (4) various flow-related phenomena. These are discussed and compared with other predictions 
reported in the literature. The adequacy or need for improvement in the spray /combustion models and the need 
for incorporating an appropriate turbulence model are also discussed. 


17. Key Words (Suggested by Author(s)) 

Rotary engine 
Combustion 
Fuel injection 

Computational fluid mechanics 

19. Security Classif. (of this report) 
Unclassified 


18. Distribution Statement 

Unclassified - Unlimited 
Subject Category 34 


20. Security Classif. (of this page) 
Unclassified 


21 . No. of pages 

24 


22. Price* 


NASA FORM 1626 OCT 86 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 






